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^ ■ Abstract 

CN . In this work we extend the JIMWLK formalism to the two-source problem. The S-matrix for the 

forward scattering can be written in a double functional integral representation which involves 

■ the usual functional integral for the gluon field and the spin path integral for the external color 
sources. Modifications needed in the light-cone gauge are discussed. Using our source term we 
compute the produced gluon field and discuss the duality of the high energy evolution kernel in 

t~^. ' the pA collision. 
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1 Introduction 



Since the pioneering work [1], impressive progress has been made in the semiclassical 
description of high energy scattering in QCD. It has been argued that, when accelerated to 
very high energy, all hadrons and nuclei attain a universal form of matter, the Color Glass 
Condensate (CGC). A hadron in the infinite momentum frame is pictured as a collection 
of random, classical charges p which produce strong color fields ~ 1/g (collective 
excitation of the small-x partons). The key issues of high energy scattering, such as 
the BFKL growth of the gluon number [2] and its eventual saturation [3], are effectively 
encoded in the quantum evolution of the CGC as summarized in the B-JIMWLK equation 
[4-6]. 

In the JIMWLK formalism, color charges p are treated as static, that is, independent of 
x + for a right-mover (the 'target'). The rationale behind this assumption is well-known, 
and this is where the name "glass" comes from: Viewed from the natural time scale of 
the small-x partons, classical charges representing high-x partons are 'frozen' due to 
Lorentz dilation. Implicit in this picture is an assumption that the effect of the left-mover 
(the 'projectile') is negligible and one is effectively dealing with a single source problem. 
Indeed, the original derivation of the JIMWLK equation [5, 6] makes no reference to the 
projectile. 
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However, in the recent surge of activity to go beyond the B-JIMWLK equation (the 
Pomeron loop physics), it has become clear that this initial picture must be modified at 
important levels. In Ref. [7] it has been recognized that in order to include a new effect 
which is missing [8] in the JIMWLK equation, one has to abandon the assumption of x + - 
independence. This effect, the gluon Bremsstrahlung, probes the multi-point correlation 
inside the target in the x + direction, and therefore it forces one to consider p as explicitly 
x + -dependent variables. Another observation, due to [9] , is that in the dilute regime where 
a hadron develops gluon number fluctuations, the non-commutativity of color charges 
should be taken into account. Again this is something that was neglected in the JIMWLK 
formalism from the outset for an obvious reason that commutators are suppressed [p, p] <C 
pp when a hadron is at saturation p ~ l/g 2 - 1 

These two apparently different issues, the x + -dependence of p and the non-commutativity, 
are nevertheless deeply related through the process of the spin path integral [10]. In a 
previous publication [11], we have fully worked out this connection in the single source 
problem. In this paper, we extend this analysis to the case of the two-source problem, 
which we propose as a model for the collision of two saturated hadrons (nuclei). In fact, 
there is already a CGC-based approach to the AA collision [12, 13] in a special gauge 
which is useful under the assumption that the sources are strictly proportional to the delta 
function in the longitudinal direction. In prospect of applications to quantum evolution, 
we relax this assumption and formulate the problem in general gauges. Our main goal 
is to set up a simple framework to study the scattering process at any stage of its high 
energy evolution, all the way from the gluon number fluctuation in the early stages 2 to 
the quantum evolution of the AA collision. 



In Section 2, we derive the source term for the hadron-hadron collision symmetrically 
for each of the hadrons starting from the eikonal formula of the QCD ^-matrix. The 
main ingredient is the path integration of color matrices. The result, Eq. (2.6), has an 
obvious problem with the light-cone gauge, and we suggest an alternative source term 
specific to this gauge in Section 3. At the end of Section 3, we compare our approach with 
related works by Lipatov [15] and by Balitsky [16]. Sections 4 and 5 are applications. In 
Section 4, we outline how to compute the produced gluon field using our formulation. A 
known result for the pA collision is reproduced. In Section 5, we discuss the duality of the 
Bremsstrahlung Hamiltonian [7, 9, 17]. In Appendix A, we give a brief review of the spin 
path integral which is used throughout this paper. 



However, even when p is large, quantum evolution of the CGC is actually sensitive to the 
commutator. See, Section 2.2 of [11]. 

2 For a recent discussion on the importance of fluctuations in the dilute regime, see [14]. 
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2 Eikonal scattering amplitude 



Let us start with high energy quark-quark scattering. In the eikonal approximation, the 
S-matrix for the forward scattering in the impact parameter space is given by a product 
of two Wilson lines [18] 

S(Y) = (V^(x)W(y)) = J [DA 11 ] V\x)W {y)e iS ^ A \ (2.1) 

where Y is the rapidity and x — y is the impact parameter. Letters in boldface denote 
two-dimensional transverse vectors. V and W are Wilson lines along the trajectories of 
colliding quarks 

V\x) = Pexp (ig J°° dx-A+(x-,x)t a ^j , (2.2) 
W(x) = Pexp (ig J°° dx + A~(x + ,x)t a ^j , (2.3) 

where P is the path-ordering operator and t a are the color matrices in the fundamental 
representation. In reality quarks have mass, so the trajectories are slightly tilted from the 
light-cone axes by a rapidity dependent angle. But for notational simplicity, we keep the 
deviation from the light-cone implicit. Also, gauge fixing (e.g., the Feynman gauge) is 
understood in Eq. (2.1) and we suppress the matrix indices of the Wilson lines. 

In QED, Eq. (2.1) can be exactly evaluated if one neglects Fermion loops [19]. In QCD, on 
the other hand, it is a very difficult problem and one needs some approximation schemes. 
One of the popular approaches is the semiclassical method where one reduces the problem 
to solving the classical equation of motion. While the effectiveness of such a method is 
not immediately obvious for the quark-quark scattering, it is arguably the most realistic 
approach for the collision of two saturated hadrons/nuclei. An obstacle to the classical 
description of Eq. (2.1) is the non-commutative color matrices t a in the Wilson lines. 
However, one can convert these matrices into c-numbers using the method of spin path 
integral [10] 

W(a;) = J [Dp(x + )] e i5 wz[p0z + )M3 / dx+ P a (x+)A-(x+,x) ^ ^.4) 

where the color charges p(x + ) are commutative, but time dependent. Swz[p] is a Wess- 
Zumino term. In Apppendix, we give a quick derivation (and a more precise form) of 
Eq. (2.4). In the context of high energy scattering, the Wess-Zumino term was first in- 
troduced in [9]. Its origin was clarified later in [11]. Using the formula Eq. (2.4) for both 
W and V, one obtains an alternative expression of the S'-matrix (the subscript L(R) is 
for the left-(right-)mover) 

S Y = J[Dp R (x + )Dp L (x-)]e iSwMx+)]+iSwz[pL(x ~ )] 

X J [DA 11 ] g i5, YM[^l]-«S / dx~ pl(x-)A+ (x~ ,y)-ig J dx+ p^(x + )Aa (x + ,x) ^ 5) 
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In Eq. (2.5), two independent functional integrals are involved; one for the gauge field and 
one for the color charges. Thanks to the latter, the external non-commutative charges are 
treated as the usual c-number source term / dxJ^A^ as in the QCD generating functional. 

We now generalize the quark-quark scattering formula Eq. (2.5) to scattering of more 
complicated objects. Following the original suggestion [1], we model a hadron or a nucleus 
at very high energy as a bunch of partons belonging to various representations of the color 
group and continuously distributed in the transverse plane as well as in the longitudinal 
(i.e., x~ for the right-mover) direction. For a thick nucleus, or a highly evolved hadron, 
the representations can be very large. Each parton propagates eikonally and leaves a 
Wilson line in the corresponding representation in its wake. Thus we propose the following 
representation of the ^-matrix 

S Y = ([lVl(x i )I[W(y j )) 

= J [Dp L (x-)Dp R (x + )]Z[p L ]Z[p R ] e i5 wz[PR(z+)]+iSwz[PL(z-)] 

X J [DA 11 ] e iS YM[AM# J d 4 xp- L A+-ig j d A xp\A~ ^ 6) 

Note that integrals in the source terms are four-dimensional. Viewed from a distance, the 
distribution of charges is localized in the longitudinal direction and may be regarded 
as a delta function in the first approximation Pr(x) ~ 5(x~)p R (x + , x) and Pl(x) ~ 
5(x + )p L (x~ , x), though it is known that regularization is needed for certain applications. 
For a 'dilute' hadron which has not yet gone through high energy evolution, p is of order 1, 
whereas for a 'dense' hadron, p can be as large as ~ 1/g 2 characteristic of saturation. Z[p] 
is a weight function which describes the probability distribution of partons at x^ — > — oo 
like the one in the dipole model [20]. 

In practice one is typically interested in the scattering of globally color neutral objects. 
In this case it is possible to recast Eq. (2.6) in a gauge invariant form. One can construct 
a closed loop out of a pair of Wilson lines corresponding to a charge and an anti-charge 
at different transverse points by attaching transverse Wilson lines U at x~ = ±oo 

V\xi)V(xj) — > U(xj,Xi,x~ = oo)V^ (xi)U(xi, Xj, x~ = —oo)V(xj). (2.7) 

The path connecting £Cj and Xj is arbitrary, and we choose it to extend to spatial infinity 
and then come back 

U(xj, Xi) = C/ f (oo, Xj)U{po, oo')f/(oo', Xi) (2.8) 

This procedure is repeated for all existing color singlet units (dipoles, quadrupoles, etc.) in 
the projectile, and similarly for the target. Then the S'-matrix Eq. (2.6) acquires a gauge 
invariant meaning. Indeed, since the field strength Fij vanishes at infinity x ± = ±oo, it 
is possible to gauge away transverse fields A i {x ± = ±oo,x). With this gauge choice a 
Wilson loop reduce to a product of light-like Wilson lines as in Eq. (2.6). However, this 
also means that in generic gauges Eq. (2.6) should be understood that at each end point 
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of Wilson lines a transverse Wilson line extending to infinity U(oo, x) is attached. 3 

Instead of using the path integral formula Eq. (2.4), one could trivially rewrite the Wilson 
lines in the first line of Eq. (2.6) as 

W(x) = J [Dp R (x + )]e- i9 I dx+ ^ x+ ^ x+ ^Pe W ^jT (fc+^A^ 5[p R (x + )], 

(2.10) 

and similarly for V\x). This is in accord with the suggestion in [21] that the weight 
function must contain Wilson lines of the functional derivative, S/Sp, to ensure non- 
commutativity. We note that the derivative Wilson line turned out to be useful for 
the generalized dipole model [22-24] (and in other context [25]). Both Eq. (2.4) and 
Eq. (2.10) are exact rewritings of a Wilson line. However, the meaning of the measure 
[Dp] is completely different in the two cases. Eq. (2.4) is the genuine path integral of the 
'spin operator' t a — > p a and therefore the measure [Dp] respects the compactness of the 
underlying group (see Appendix). One can speak of its saddle point, gauge transforma- 
tion, etc. [The saddle point approximation will be important in Section 4.] On the other 
hand, p R (x + ) in Eq. (2.10) is a dummy variable which comes from the insertion of unity 
l = f[Dp] e-^ A S[p]. 



3 Light cone gauge 



For the description of high energy scattering, it is often very convenient to work in the 
light-cone gauges, A + = or A~ = 0. The source terms in these gauges assume a peculiar 
form and deserve a separate study. For definiteness, we choose A + = 0. In this gauge, 
= 1, and as one can see from Eq. (2.7) the transverse Wilson lines represent the 
dominant interaction with the left-mover. 

Instead of setting A + = in Eq. (2.7) from the beginning, let us reach this gauge via 
gauge rotation from other gauges in which Ai(x~ = ±oo) = 0, that is, Eq. (2.6) is valid 
without insertions of transverse Wilson lines. Such a detour is useful since the relevant 
gauge rotation often appears in the literature [6,7]. Besides, it prescribes the boundary 

3 Incidentally, we note that for a closed loop it is possible to rewrite the Wess-Zumino term in 
a more elegant way than that given in the appendix. In the case of color SU(2), 

Pexp (^g j> dxPA*?^ = j [Dp(t, s)] exp (^-ig j> dt^A^p a (t) + iJ j> dt J ds(d t p x d s p) ■ jjj 

(2.9) 

where t is the coordinate along the contour (the 'ordering variable') and s is the radial coordinate 
of the disc spanned by the contour. p(t, s) is a continuation of pit) inside the disc such that at 
the boundary, p(t,s = 1) = p(t). This (1 + 1) -dimensional Wess-Zumino term (often called the 
Polyakov spin factor) is the form suggested in [9]. 
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condition of the transverse field at x — > ±00. The A + = gauge can be achieved by the 
following gauge rotation 

A^ 1 — > A^ 1 = VA^V ] + -Vd»V\ (3.1) 

9 

where (c.f., Eq. (2.2)) 

V\x) = Pexp (ig J X dz-A+{x + ,z-,x)t a ^J . (3.2) 

Then the time-like Wilson line in the original gauge can be rewritten as a transverse 
Wilson line in the A + = gauge 

V\x) = Pexp [ig J°°dz l A^{x- = 00, z)t a ^j = e" 1 ^*", (3.3) 

where the integration path is arbitrary since Ai is pure gauge. To make contact with the 
light-cone gauge perturbation theory we note that to leading order in g, Ai ^ <9j0, and 
Eq. (3.3) can be written as 

V\x) » exp (igt a ^A$(x- = oo^J . (3.4) 

The coupling ~ (d' l /'V 2 )Ai is precisely the structure one finds from the direct perturbative 
analysis in the light-cone gauge (see, e.g., Sect. 3.2 of [27]). Eq. (3.4), or its continuum 
generalization t a — > — / dx + dxp a L (p being non-commutative) is the source term that is 
useful for our purposes. Due to the long range nature of the light-cone gauge propagator 
~ l/p + , sources at infinity x~ = 00 cannot be neglected. The scattering in this gauge is 
pictured as a left-moving quark changing directions at x~ = 00 to the transverse direction, 
and the interaction with this quark entirely comes from this part of the propagation. 

The usual assumption of saturation amounts to postulating that pL ~ I/*? 2 and dropping 
non-commutativity by hand Pl — > Pl- 4 In place of Eq. (2.6), we now have (omitting the 
tilde on A^) 

J [Dp R (x + )D PL ]Z[p R ]Z[p L ] e ^wz[p R (*+)] 



X 



j A+= [DA»] exp [iS YM [A] -ig J d i xp K A--ig J^ dx +dxp a L —AA . (3.5) 



The role of the source term ~ p L is to generate a classical field of the left-mover (see, [16] 
for a construction of source terms from this perspective) 

A d = -^p L (x + ,x). (3.6) 
4 Note that the loss of non-commutativity also justifies the approximation di4> ~ A; L . 
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Indeed, it precisely cancels the boundary term of the Yang-Mills action when expanded 
around Eq. (3.6) 



In general, by choosing a prescription for the singularity l/p + of the light-cone gauge 
propagator, or equivalently, by fixing the residual gauge freedom of the A + = gauge, one 
can require either A l (x~ = oo) = or A l (x~ = — oo) = [6,26,27]. Eq. (3.5) corresponds 
to the latter (the final state interaction) in accordance with the choice Eq. (3.2). Of course, 
one could have chosen the other boundary condition Ai(x~ = oo) = (the initial state 
interaction) by changing — oo — > oo in Eq. (3.2). 

The above procedure is valid only when the effect of the right-mover is negligible, pr ~ 1. 
If the right-mover is at saturation, p^ ~ 1/g 2 , if would create a strong classical back- 
ground field ~ 1/g and the source term must be modified accordingly. We assume that, 
in the original gauge, 



where a ~ 1/g is the classical field of the right-mover. We then transform to the A + = 
gauge in two steps. First we rotate away the classical component 




))• 



(3.7) 



A + = a + 5A + , 



(3.8) 



A" -> A" = VqA^vJ + -V d»Vj, 



(3.9) 
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where 




(3.10) 



In this gauge, 



A+ = V 5A + Vj, 



A\x- = oo) = -Vq&vJ = B\ 



(3.11) 



Then we remove the remaining weak field A + by the gauge rotation with 




(3.12) 
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The transverse field at infinity becomes 

Ai(x- = 00)^ 5VBi5V ] + l -5Vdi5V ] ^B { + Drf, (3.13) 

where A = di — igBi. Assuming that the left-mover is also saturated, we can let t£ — > 
— / dx + dxpf j and express in the original gauge in terms of the boundary transverse 
field in the light-cone gauge Eq. (3.13). 

V* = V^5V ] = e ~ i9X e~ i9<t> — > e ig J ^ e ig J p £<t> a 

= e -i j DiA-DiX-i J D i A^ l D i <f> 

« exp (i j dx + dxF~ l i {B i + Drf)) , (3.14) 

where is the field strength constructed from the two classical background fields 
(Bi, A d ) with D V F"{~ = D i D t A~ l = gp L (c.f., Eq. (3.6)). In the last line, we have used 
diX ~ Bi and neglected the commutator [-B«,x] ~ [Bi(x), J x dz l Bi(y)] assuming that Bi 
does not change much over the typical separation of a charge and an anti-charge in a 
dense projectile. Eq. (3.14) is the desired source term for the left-mover with the field 
of the right-mover taken into account. We see from Eq. (3.7) that its linear part pre- 
cisely cancels the boundary term of the Yang-Mills action when expanded around the 
background field A» -> B i + + A». (In fact, it also cancels S Y m[A^}. See, Eq. (5.6) 
below.) 



Summarizing, Eq. (2.6) and Eqs. (3.5), (3.14) are our setup for the two-source problem. In 
the semiclassical framework one first solves the Yang-Mills equation with sources obtained 
by varying with A p in Eq. (2.6) 

D V F V » = 5^p R (x + ) + 5»-p L (x-), D- PR + D + PL = 0. (3.15) 

The second equation, the current conservation law, is a consequence of the first. Although 
Eq. (3.15) is frequently written down in the literature of the CGC, the precise meaning of 
the time dependence of charges, or from which source term the equation has been derived 
is not always articulated. This latter point deserves further comments. Assuming that the 
currents are separately conserved, 

D+ PL = 0, D dPR = 0, (3.16) 

in our approach we can integrate these equations with an arbitrary boundary condition, 
in particular, the retarded boundary condition 

Pl(x-) = V^x-)p L (-oo)V(x~), 
Pr(x + ) = W(x + )p R (-oo)W\x + ), (3.17) 
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which is usually employed for the gluon production problem. |V(x - ) (W(x + )) is the same 
as in Eq. (2.2) (Eq. (2.3)) except that the upper limit of the integration is x~ (:r + ).] On 
the other hand, currents derived from effective source terms of the JIMWLK formalism 
such as Sw = / d 3 xtr[p R W] and S\y = / d 3 xtr[p R In W\ (W being the adjoint Wilson line) 
[5,6,28,29] via J + = 5Sw/5A~ are covariantly conserved, D~J + = 0, but one cannot 
freely choose the boundary condition to this differential equation because it is already 
prescribed in the structure of S\y- I n fact, it has been pointed out in [6] that any effective 
source term of the form Sy/[p, W] is incompatible with the retarded boundary condition. 
Our source term with the Wess-Zumino term is not subject to such a constraint. 

In practice, one has to resort to some sort of perturbative expansion to solve Eq. (3.15). 
The Yang-Mills action plus source terms evaluated at this saddle point is often referred 
to as the effective action. Clearly, the choice of the source term and how one expands the 
saddle point solution are the key issues. We shall work out the case of the pA collision 
in Section 5 and leave the analysis of the AA collision to future works. [This was partly 
done in [7].] In the reminder of this section we draw a comparison (albeit superficially) 
between our approach and two previous approaches in the literature. Our main concern 
here is the form of the source term. 

(i) In the context of quark-quark scattering, Lipatov [15] constructed an effective action 
describing the interaction between the s-channel 'physical' gluons V 1 and the t-channel 
'reggeized' gluons A ± in a given rapidity interval 

S cS [v] = S Y m[v] -—I dx-dxtT[W[v-]V 2 A + ] - — f dx + dxtr[V^v + \V 2 A-\ + ■■■ , 
9 J 9 J 

(3.18) 

where v^{x) = A + (x~, x) + A~(x + , x) + V fl (x) is the total gauge field. Eq. (3.18) describes 
quasi-elastic processes (linear terms in A^), the production of gluon clusters in the multi- 
regge kinematics (quadratic terms in A^), and so on. One can integrate out the V^ 1 field by 
solving the classical equation of motion 5S eS /5v^ = perturbatively in powers of g. The 
effective action evaluated at this solution S , e fr[vci[^4 ± ]] contains all possible reggeon number 
changing vertices in the t-channel (A + ) m — > (A~) n in a given rapidity interval. Clearly, 
there is a close correspondence with the effective action for the two-source problem in the 
CGC approach. One notices that the structure 5 of the interaction terms in Eq. (3.18) is 
identical to the source term of the JIMWLK formalism 5 W oc / d 3 xtr[p R W] [5,6]. It is 
thus tempting to identify the reggeon fields in [15] with the classical color charges of [5, 6] 

V 2 A + <-> gp R (x~,x), 

V 2 A- <-> g PL (x + ,x). (3.19) 

Eq. (3.19) may be related to the technicality in [15] that the reggeon fields are invariant 
under local gauge transformations. It would be important to study how the effective action 

5 Wilson lines in Eq. (3.18) arise after determining the boundary condition for the inverse 
operator c^ 1 in the multi-gluon production vertices. A more fundamental expression of the 
interaction term can be found in [15]. 
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in the Wilson line representation [7, 16] is matched with the perturabative expansion in 
[15]. 



(ii) Instead of introducing color charges, Balitsky described the colliding hadrons as shock 
waves, i.e., fixed gauge field configurations [16]. Eq. (2.6) may be compared with the twice 
rapidity-factorized formula for the scattering amplitude 

iT Y (p a:Pb - p' a ,p' b ) = J [DB"DC"\ e iS ^ +iS ^J(p a )J(p b )J(p' a )J(p' b ) 

x J [DA"} e iSYM[A]+i I fofiWW+t I (3.20) 

where J's are the currents of external hadrons and 

/oo 
dx + W\x + )F +i (x + )W(x + ), 
-00 

/oo 
dx-V{x-)F^{x-)V\x-). (3.21) 
-00 

The source terms in Eq. (3.20) perfectly makes sense in any gauge including the light-cone 
gauges. The gauge fields A, B, C are ordered in rapidity, Yb > Y\ > Ya > Y 2 > Yc and 
the Wilson lines in / and g are actually not strictly light-like but along the direction of 
the rapidity divider Y\ and Y 2 , respectively. Viewed from the A" field, Eq. (3.20) is a two- 
source problem just like Eq. (2.6). It has been suggested [16] that the classical Yang-Mills 
equation with sources ~ Sfi/SA~, Sgi/SA + can be formally solved via the commutator 
expansion. 



4 Gluon production 



Though Eqs. (2.6), (3.5) were derived as the ^-matrix for the elastic scattering, proper 
specification of the source terms should allow one to compute the particle production in 
hadron collisions as already implied around Eq. (3.17). In this section we demonstrate 
that this is indeed the case. The produced gluon field (the inclusive gluon production) 
can be computed via straightforward manipulations of the functional average of the gauge 
field in the presence of sources (Hereafter we often omit '/ <i 4 x' to simplify the appearance 
of equations. ) 

(A"(x)) = y*[ J Dp L (a;-) J Dp R (a; + )]e i5wz[m(:E+)]+i5wz[pL( ^ )1 

x J[DA"]e lSYM[A] - ipnA ~- l ^ A+ A"(x), (4.1) 

with the premise that the last propagator contracted with A"(x) should be the retarded 
one. Typically one is interested in the behavior in the forward light-cone or* 1 > 0. To 
see how Eq. (4.1) works in practice, below we consider two examples, the pp and the pA 
collisions. In both cases we shall illustrate the nontrivial roles played by the Wess-Zumino 
term. 
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4-1 The pp collision 



When both charges are small p R / L ~ 1 (the pp collision), (A fl (x)) can be computed 
in perturbation theory by expanding g dependent terms (source terms as well as the 
interaction terms in the Yang-Mills action). One of the contributions to g 3 is given by 

t^l{ pR A-( w )p R A-(y)p L A + (z)A' i (x)) = Hg) 3 G^ (x,w)G ^y, z)(p R (w)p b R (y) P Uz)), 

(4.2) 

where Gq V is the free gluon propagator. The product of charges should be handled a little 
carefully. Due to the presence of the Wess-Zumino term, the following identity holds (c.f., 
Eq. (A.14)) 

(p R (w + , w)p R (y + , y)> = 6{w+ - y + )p a R (w)p b R (y) + 6{y+ - w + )p b R (y)p R (w), (4.3) 
where p R is the non-commutative, but x + -independent charge of the right-moving proton 
«(y) " Pr^PrM = -if^fikivMw - y). (4.4) 

(p L and p R commute.) We decompose the product in Eq. (4.3) into symmetric and asym- 
metric parts to get 

^f ab \e{w + - y + ) - %+ - w + ))p R (y)S(w - y) + \{p R {w), p b R (y)}. (4.5) 

The symmetric part does not contribute when we put the detected gluon on-shell, so it 
can be discarded. Other terms are similarly computed and the result is identical to [30]. 6 

4-2 The pA collision 

Next we consider the pA collision with the proton moving to the right. By this we mean 
that p R ~ 1 and p L ~ l/g 2 - This problem has been extensively studied in the literature [26, 
31-33], so we merely intend to show how things work in our formulation. It is convenient 
to choose the A + = gauge. The produced gluon field is (see Eq.(3.5)) 

(A' t (x))= [ [Dp R (x + )} I [DA"} e iS ^[A]-i 9 pnA-+i9^ Pl ,A i ( x +,o , X )+iS wz [ P n(x+)] A ^,y 

Ja+=o J 

(4.6) 

Unlike in the pp collision, one cannot expand in powers of p-^A 1 since pl is large. However, 
this term can be removed (see Eq. (3.7)) by making the following shift 

A" = ^ + a", (4.7) 

6 Except for the e-prescription of internal Fermion propagators 6{w + — y + )^> ^(9(w + — y + ) — 
6{y + — w + )) which, however, does not affect the result if the outgoing gluon is put on-shell. 
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with A cl given in Eq. (3.6). Then Eq. (4.6) becomes 

(A»(x)) = [ \Dp R {x + )} e ^wz[PR(-+)] f [ Da ^ e i^D-i[A-]a"-i 9 pn(A-+a-)^ A - + ^ 
J Ja+=0 

(4.8) 

Let us focus on the transverse field (A l (x)), since this is the only component which con- 
tribute to the gluon production in this gauge. Since p R is small, we expand only to first 
order in p R a~ . 

(A\x)) = j[Dp R (x + )} e < Swz[p R (x+)] | [jDa M] e i^D-^A-^-i 9PK A-^_ ig ^ J pKari<y y^ 

= -ig [ [Dp R {x+)]e^^ + ^^ A * ( dy + dyG*-(x,y)[A d ]p R (y+,y), 
J Jy~=0 

(4.9) 

where G^LA^] is the gluon propagator in the a + = gauge in the presence of the back- 
ground field A~ Y . In Eq. (4.9), we can safely set p R (y + ,y~,y) = p R {y + ,y)5(y~). It is 
possible to do the saddle point approximation for the path integral Dp R (x + ) since a large 
factor A~\ multiplies p R . The saddle point equation is nothing but the current conservation 
law [11] 

D- lPR (x + ) = 0, 

- Pn(y + ,y) = W(y + ,y)p R (y) = Pexp [ig J^dz + A^z + , y)^j p R (y). (4.10) 

At the saddle point, Swz[pr(^ + )] — 9PrA~ x = [11]. The background field propagator can 
be extracted from [1, 6] 

y) = J dp + [Go(x-y,p + )9(y + )d^ 

+ 2ip + J dzG {x-z,p + )5(z + )G Q {z-y,p + )6(-y + )w)jdl l 
= 2nj dz-6(z-)G (x -fr x -- y -- Z -)9(y + )d* 

+ 4tt J dzG (x - z)5(z-)G (z - y)6{-y + )W z d^ (4.11) 

where x = (x + , x) and 

W y = Pexp (ig dy + A A {y + , </)) . (4.12) 

In Eq. (4.11), we have dropped terms which do not contribute at x + > 0. means that 
the derivative acts on functions on its left. Go is the retarded free scalar propagator. In 
the second equality, we have chosen the light-cone gauge prescription l/p + — > 1/G» + +ie) 
to go with the boundary condition we have chosen. After the substitution of Eqs. (4.10) 
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and (4.11), Eq. (4.9) takes the form 

(A\x)) = -2mg J ^ dy + dy (| d Z -9{ Z -)G {x -y,x~- z-)9(y + )d^ 

+ 2 J dzG (x-z)5(z + )G (z-y)6(-y + )W z dl^W(y + ,y)p R (y) 
= 2mg f d 4 y9(y~)9(y-)G (x-y)dl(W y p R (y)) 



+ 



*mg J_ =Q d y +d y J dzG (x-z)5(z + )G (z-y)9(-y+)W z dlp R (y), 



(4.13) 



where we used W{y + ) = 9(y + )W + 9(—y + ) since the support of p^ (hence that of A cl ) is 
small in the x + direction. One can simplify the second term in (4.13) 

L, _,/f + * G oO - V)0(-V*)W z d' y f >R (v) 

J (27r) 2p p + — p + ip + ep — le y 

r dr> + dn e~ iv ^ ~ z ~ ~+ i p(. z -y) 

V 

and the result for the produced gluon field is 



= S(z-)W z ^p R (z) (4.14) 



(A\x))=2mg[ / d*y9(y- )9(y + )G (x - y)d> v (W vPR (v)) 

di 



+ 2 J dzG (x - z)5{z + )5{z-)W z ^- 2 p R {z)y (4.15) 
in agreement with [33]. 



5 Revisiting the Bremsstrahlung Hamiltonian 



In [7], the dual partner of the JIMWLK Hamiltonian-the Bremsstrahlung Hamiltonian 
-f^BREM-was derived. [See also an earlier derivation in [9] .] However, its concrete usage was 
not clearly explained there, nor was the source term specified. In this section, we complete 
the calculation of [7] by using the source terms obtained in Section 3. As we shall see, this 
step is crucial for the proper discussion of duality [7, 17] of the evolution kernel. We also 
show how our formulae of the ^-matrix containing two eikonal phases (e.g., Eq. (3.5)) 
reduce to the CGC-type formula (see Eq. (5.11) below) of the S'-matrix containing only 
one eikonal phase. 

We consider exactly the same situation as in Section 4.2, namely, scattering of a dilute 
right-mover (p R ~ 1) and a dense left-mover (p L ~ V^ 2 ) i n the A + = gauge. The 
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forward S'-matrix is (c.f., Eqs. (3.5), (3.14)) 

S Y = J [Dp R (x + )Dp L ] J [/}^] e ^YM[A]- lWR (a ; +)A-+iF c 7M l (x-=oo)+^wz[PR(^ + )] > ( 51 ) 

We expand A p around a sum of solutions to the single source problem 

A^ = B i (x + , x~, x) + A~(x + , x), (5.2) 

where 

-D i d + B i = gp R (x + ), -V 2 A-=gp L . (5.3) 

Note that = 0. In [7], the current conservation law D~ x p R = was imposed. However, 
the first equation of Eq. (5.3) has a solution of the pure gauge type F C Y = irrespective of 
the conservation law, and below we do not impose it. After a shift A^ = A^ + a M , Eq. (5.1) 
reads 

S Y ~ J [Dp R (x + )Dp h ] e l5 wz[p R (x+)] 

x J [Da"] exp (-i 9PR A- + iD^a* + Vg^K + •••)■ (5.4) 

Notice that the leading term is given by a single eikonal phase e~ pA due to nontrivial 
cancellation between 

S YM [A d ] = F c fF^ » —d + B i D i A~ = -D,d + B l A~ = gp R A~, (5.5) 

and 

F d i B i (x + , oo, x) « -D i A~d + B i = A~D i d + B i = -gp R A~, (5.6) 

where we used B l (x~) = 6(x~)B l (oo). Since \B l \ <C \A~\, we have kept only terms linear 
in B\ It is crucial to specify the correct source terms to see this cancellation. 

Performing Gaussian integration in Eq. (5.4), and reinstating the weight functions, we 
obtain 

S Y = J [Dp R (x + )Dp L }Z Y _ Yo [p R }Z Yo [p L ] 

x exp (iS wz [p R (x + )] - igp R A- + i5Y S BREM [p R (±oo) , W]) , (5.7) 

where SY is the rapidity interval carried by the a" field. As indicated, the right-mover 
carries rapidity Y — Y and the left-mover, Y . Sbrem is an effective action describing 
gluon Bremsstrahlung whose explicit form can be found in [7]. It contains charges p R at 
x + = ±oo and the adjoint Wilson line W. From Eq. (5.7) one can deduce the evolution 
equation 

[Dp R (x + )Dp L ]Z Y _ Yo [p R ]Z Yo [p L ] iS BREM [p R (±oo), W\ exp (iS wz [p R (x + )] - igp R A~ 

(5.8) 



dY 
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There are two ways to proceed from Eq. (5.8): 

(i) The eikonal factor e~ %9PRA allows one to replace pr(±oo) — > i5/5A~(±oc) in Sbrem 
to convert the action into an operator, the JIMWLK Hamiltonian 



25brem[p±oo, e 



-H 



JIMWLK 



5 



5A 



±00 



(5.9) 



After this is done, one can use the path integral formula Eq. (2.4) in the backward direction 
and get 



dS Y 



= - [ Dp L Z Yo [p L }H J1MWLK ^^,W(A-(x + )) (UW( Xi )) Y _ Yo , (5.10) 



with 



S Y = f Dp L Z Yi) [p L ] (n^)>y-y 0l 
J i 

(IJ^i)) = / Dp R Z[p R ]l[W(x^. 



(5.11) 



The differential operator ff jimwlk can act on either (left or right) direction through the 
integration by parts. When it acts on Z[pi], Eq. (5.10) describes the evolution of the 
left-moving Color Glass Condensate, namely, the JIMWLK equation. When it acts on 
the right, Eq. (5.10) describes the evolution of Wilson lines, the Balitsky equation. In the 
latter case it may be more natural to put the energy (Y — Y ) dependence in the slope of 
the Wilson lines as originally done in [4]. 

(ii) Returning to Eq. (5.8) this time we replace A~(x + ) — > i5/5p R (x + ) in the adjoint 
Wilson line W in Sbrem to obtain the Bremsstrahlung Hamiltonian 



^Sbrem[p±<x>, e 



iA- 



6_ 

-#BREM[P±oo,e Sp\. 



(5.12) 



After this is done, -Hbrem has no reference to the left-mover anymore, so we can drop the 
averaging over Z[pi] 



_d_ 

W 



Y[W(xi)) = -J Dp R (x + )Z[p R ]e lS ^ 



(*+)] 



x ff , 



BREM 



p R (±oo), w(j-^)\ exp (-igp R A-) . 



(5.13) 



Functionally differentiating both sides with respect to A (x + ) arbitrary times and setting 
A~ = 0, we obtain an equation for the correlator of non-commutative charges 



_d_ 

W 



<PrVrVr c • • • ) 



x ff , 



BREM 



Dp R {x + )e lS ^ p ^Z[p R ] 
5 



PR 



±00), w(— — — ) 
{ -5p R (x+)^ 



(5.14) 
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where the 'ordering variables' xj satisfies the constraint x± > x\ > x% > ■ ■ ■ but 
otherwise arbitrary. Eq. (5.14) coincides with the equations derived in [9]. Our derivation 
unambiguously shows that the correct ordering variable is x + , and not x~ as claimed in 
[9,21]. 

The above calculations nicely exemplify the duality property of the evolution kernel [7, 17]. 
The effective action <Sbrem can simultaneously describe the gluon saturation of the left- 
mover (Eq. (5.10)) and the gluon Bremsstrahlung of the right-mover (Eq. (5.14)) in 
asymmetric scattering. 



6 Conclusion 

In this paper we considered scattering of two saturated hadrons in the functional integral 
framework. We have constructed the source term for each of the hadrons including the 
case of the light-cone gauge. Our formulation is applicable to colliding hadrons at any 
stage of their high energy evolution, from the dilute regime p ~ 1 to the dense regime 
p ~ l/g 2 - Once one knows the proper source term, issues such as gluon production and 
the evolution of the S'-matrix become more transparent and straightforward to deal with 
as demonstrated for the pA collision. We hope to address the problem of the AA collisions 
in future work. 
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A Path integral of a spin 

In this appendix we give a very quick review of the standard path integral procedure for 
a spin operator. Rigorous treatment can be found in the literature [10]. For the sake of 
simplicity, we confine ourselves to the SU(2) spin. 

For an ordinary quantum mechanical system endowed with a Hamiltonian H [q, p] (q is the 
canonical coordinate and p is its conjugate momentum), the path integral representation 
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of the transition amplitude is elementary, and is given by 

(q f \e-* TH ^\q t )= [Dq{t)Dp{t)\ exp [i £ dt(p(t)q(t) - H\p(t), g(t)])) . (A.l) 

If the Hamiltonian contains spin operators, the path integral becomes less straightforward. 
As the simplest example, consider a non-relativistic particle with the spin operators S a 
(a = x,y,z) in the J-representation (J = 1/2, 1,3/2, • • •) of SU(2) placed in a constant 

— * 

magnetic field B pointing to the z-direction. The interaction Hamiltonian is 

H = -S ■ B = -S Z B. (A.2) 

The spin operators satisfy the usual commutation relation [S a , Sb] = ie a b c S c . Our goal is 
to write down a path integral formula similar to Eq. (A.l) for the Hamiltonian Eq. (A.2). 
We must begin with finding out what are q and p in this case. 

In the path integral representation, a quantum operator q has its classical counterpart 
q(t). To the spin operator corresponds a rotating vector on a sphere 

— * 

S{t) = Jn{t) = J(sin6> cos0, sin 9 sin 0, cos#), (A. 3) 

— * 

where n is a unit vector. Though S has three components, there are only two independent 
variables, 9 and <p. In fact, plays the role of the canonical coordinate, and the conjugate 
momentum is cos# 

q(t) = <p{t), pit) = J cos 9(t). (A.4) 
Indeed, if one imposes the Poisson bracket 

{0, Jcos#} R b. = 1, {M}p.b. = ~T~ r ~Q-> ( A - 5 ) 

J sin 9 

from Eq. (A. 3) one can easily derive the expected Poisson brackets among the components 

— * _ 

of the vector S ' 

{S a , S b }p. B . = e abc S c . (A.6) 



With the canonical coordinate and momentum at hand, we naively insert them into 
Eq. (A.l) and obtain a path integral formula for the Hamiltonian Eq. (A.2) 

r <P(T)=<j> f ( r-T , . . \ 

(6Ae- lTH{s \<p t ) = / [D(f)(t)D cos 9(t)} exp [i dt( J cos 9<f) + BJ cos 9) . (A.7) 

It turns out that Eq. (A.7) is (more or less) the correct formula, though the actual proof 
[10] requires careful manipulation of integrals. One starts with the discretized form of the 

7 As usual, this differs from the quantum commutation relation [S a , S b ] = ie abc S c by a factor i. 
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right hand side of Eq. (A. 7) 

OO -y 

Elim - — - 
n->oo (2ir) 

n=—oo V "/ 



I „A' N N-l 



£ lim —— f f f[Dcos6 k l[ Dfaexp *E (t(jcos6 k + 1 ) (pk + egJcosfl, 

(A.8) 



fc=i fc=i 



fc=i 



where T = Ne, 4>'f = 0/ + 27m = 07v, 0« = 0o> an d shows that it reduces to the left hand 
side 

e-* TH %,> = — V e i ^-^ +BT )^. (A.9) 
2ir S z =-J 

A constant 7 is needed in Eq. (A.8) for a technical reason. For half-integer J, one can 
choose 1 = \- 

The kinetic term / pq ~ / 00S6 1 • in Eq. (A. 7) is what we call 'the Wess-Zumino term'. 
To obtain the formula for a Wilson line, Eq. (2.4) in the main text, one just needs to 
make the following replacements in Eq. (A. 7) 

t^x + , B a ^gA-(x+), S a ^t a , S a (t) - -p a (:r + ), (A.10) 

and use a sloppy notation for the measure D<p{t)D cos #(£) — >■ Dp(x + ). In practice, one is 
interested in the transition amplitude between states with definite color indices Wij(x), 
rather than definite values of the conjugate variable 0. In this case one has to project onto 



these color states by integrating over the initial fa and final 0/ angles: J 27r dfaje 
The Wess-Zumino term has a simple representation in terms of a matrix 

S WZ = J J dx + cos 6<j) = U J dx + tr [r 3 Md+M% (A. 11) 

where r 3 is the third Pauli matrix and 



M = 



( 6 i& -9 -ii- \ 

cos|e*2 5 sin |e 2 4 
- sin |e 2 ; cos |e 2 



(A.12) 



Furthermore, M is related to p via 



p a = - -tr [r 3 Mr a M t ]. (A.13) 
2 

Eqs. (A. 11) and (A.13) are the forms used in [11] (there the letter S was used instead 
of M). Since Eq. (A. 7) (or Eq. (2.4)) is the genuine path integral for the spin, the non- 
commutativity of S a (or t a ) is guaranteed to be preserved. Yet, in the path integral 
representation, of course one can treat S a (t) (or p(x + )) as c-numbers. For instance, the 
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following identity holds 

(f\ t a H a H<* ■ --t^i) = J \Dp(x+)} e l5wz[ ^ +)] p ai «)p a2 «)p a3 «) • " -P a "«), 

(A.14) 

where x + 's are constrained such that x+ > x+ > x+ > ■ ■ ■ > x+ n , but otherwise arbitrary 
(only the ordering matters). If one changes the order of matrices, say, t a2 and t" 3 , one 
correctly obtains the commutator term if a ^ a 'i c p c on the right hand side. 
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